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1-1 ■ An efficient numerical algorithm to evaluate one-loop amplitudes using tensor 

! integrals is presented. In particular, it is shown by explicit calculations that for or- 

t> ' dered QCD amplitudes with a number of external legs up to 10, its performance is 

^ ■ competitive with other methods. 

O : 

in ■ 
o ■ 

on ' 1 Introduction 

o ■ 

In order to deal with the data from the experiments at LHC for the study of elementary parti- 
cles, signals and potential backgrounds for new physics have to be under control at sufficient 
accuracy [1]. In particular, hard processes with high multiplicities, involving many particles or 
partons, cannot be neglected. On top of that, such processes have to be dealt with at the next- 
to-leading order (NLO) level to, for example, reduce the scale dependence of observables and to 
have a better description of the shape of their distributions. 

An important part of a NLO calculation concerns the one-loop amplitude. Recently, im- 
pressive results have been published for one-loop QCD amplitudes for very high numbers of 
partons [2-4]. They were obtained with the so-called unitarity-approach. Originally restricted 
to analytical calculations [5-9], the potential of this method in a numerical approach became, 
after the crucial input from [10], clear with the work of [11-13] and [14]. It is considered an 
alternative to the "traditional" approach involving tensor integrals. Both approaches expand the 
one-loop amplitude in terms of a basis set of one-loop functions. In the unitarity-approach, this 
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set consists of scalar-integrals up to 4-point or 5-point functions, and it aims at determining the 
coefficients directly. In the "tensor-approach", the basis set is larger and consists of tensor inte- 
grals or their coefficients functions when expanded in terms of Lorentz-covariant objects [15-35]. 
Also these basis-functions are eventually calculated by expressing them in terms of a smaller set 
of scalar-integrals, but this happens in a, for the particular method, universal way, independent 
of the amplitude. 

Multiplicities with up to 20 partons as achieved in [2, 3] are unattainable in the tensor- 
approach because of the asymptotic computational complexity of the latter. It arises because 
the basis set contains n-point functions where n goes up to the total number of external legs of 
the amplitude. Let us make a crude comparison between the unitarity-approach and the tensor- 
approach of [30] in which the basis set consists of "normal" tensor integrals carrying explicit 
Lorentz-indices. A first step in the analysis of the computational complexity of the two methods 
is the determination of the number of coefficients to be evaluated in case ordered amplitudes 
have to be calculated. For the unitarity-approach which determines coefficients up to 4-point 
functions, it is given by 

u\ /u\ fn\ /n\ 14n+11rL 2 -2n 3 + u 4 
4) + U + (2) + (,)= 24 ' (1) 

The number of tensor integrals is, using the fact that only symmetric tensors have to be consid- 
ered of a rank not higher than the multiplicity, 

Y- Y- ( U ) C + 3 ) = ~ 384 + 592n + 203n2 + 26n3 + n4 2 n (2) 
k=i i=o \ I / 384 

Here, all tensors up to the maximal ranks have been included. The number is obviously much 
larger than the number of scalar functions. The asymptotic behavior of 2 n for the tensor integrals 
is particularly disastrous. It is a result of the expansion in terms of n-point integrals. The 
accompanying n 4 -behavior stems from the symmetric tensor components. Of course this is not 
the whole story. Also the operations to be performed in order to determine the coefficients have 
to be taken into account. For the unitarity-approach as presented in [2] for example, this leads to a 
final computational complexity of 0{n 9 ) . In this write-up, we will see that, by using recursion on 
both the tensor integrals and their coefficients, the complexity as given in Eq.© does not change. 
Although the asymptotic complexity is exponential as opposed to the polynomial complexity of 
the unitarity-approach, it may be competitive for moderate values of n, and in fact, we will see 
that in practice it is up to n = 1 0. 

Besides the computational complexity, also the numerical stability is an important issue con- 
cerning one-loop calculations. It is typically related to inverse Gram-determinants approaching 
zero. When using tensor integrals, this issue is isolated to the calculation of the tensor integrals 
themselves. In that sense, it allows for universal solutions, and several methods to achieve this 
exist. As a last resort, the calculation of the tensor integrals can be performed at higher precision 
level for phase-space points at which numerical instabilities occur, and the decision to do so can 
be made at relatively low cost. In the unitarity-approach, numerical instabilities can show up in 
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the computation of the coefficients for the scalar functions. At the moment, there is no better 
cure known than to increase the precision level for the full calculation for phase-space points at 
which numerical instabilities occur. Fortunately, also in this case the decision is relatively cheap. 

A final issue worth mentioning is the potential towards automation of the method. The 
unitarity-approach proofs to be successful in this respect because it optimally allows for the use of 
existing tree-level machinery related to the calculation of off-shell currents or "sub-amplitudes". 
We will see that the method presented here allows for the same. In particular, it is completely 
numerical and no computer algebra is involved. 

The outline of the paper is as follows. In Section [2] the tensor reduction is addressed, and 
in Section [3] tensor symmetrization, which is crucial for the efficiency of the presented algo- 
rithm. Recursive relations for one-loop amplitudes are presented in Section HI and in Section [5] 
results can be found obtained with the help of an explicit implementation of the algorithms in the 
foregoing sections. The conclusions in Section [6] finally close the paper. 

2 Tensor reduction 

Tensor integrals are usually calculated using recursive equations relating high-multiplicity and 
high-rank tensor integrals to lower-multiplicity and lower-rank ones. Tensor reduction is this 
formal process, in practice the opposite process, tensor building, is performed. The multiplicity 
n and and rank r are defined with the formula 



As the formula suggests, we consider tensor integrals defined in D dimensions, but only with 
4-dimensional components of the integration momentum in the numerator. This will lead to a 
calculation of the one-loop amplitude within the scheme of [37], which asks for a finite coun- 
terterm in order to arrive at gauge-invariant results and to cast the result into other schemes like 
't Hooft-Veltman or FDH. This finite counterterm is exactly given by the so-called C^2-term, 
showing up explicitly in the OPP unitarity-approach as part of the rational terms [13], and which 
is shown how to be determined in [38]. 

The asymptotic computational complexity given in Eq.© does not increase when the oper- 
ations needed to calculate the tensor integrals are taken into account, because for high-n each 
integral can be obtained using a fixed number of lower integrals, independent of n or r. This can 
be easily understood as follows. Using the fact that we can write 



2(Pj-Pn) • q = [(q+Pj) 2 -mf] - [(q+p n ) 2 -m n ] + m 2 - p 2 - m n + p n , (4) 



-V1 "V2 — "Vt 

n,r 



[• d D q qj 1 qp- qg 

. m D/2 nr=i[(q + Pj) 2 -™; 




(3) 



we have 



2( Pj -Pnk^; 



n-1 ,r-1 



n-1 ,r-1 



n 



+ (mf - p 2 - m 2 + p 2 ) 7\ 



.-V1-V2— -Vt-1 

n,T— 1 



(5) 
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where T^ 2 r lY 1 0) * s obtained from O^ 1 ^ "" Vr_1 by removing the j-th denominator. Choosing 
4 different vectors Pj appearing in the denominators, we get 4 relations, enough to determine 
the 4 integrals T^V 2 " with the first r — 1 Lorentz indices fixed. So 4 tensor integrals can 
be determined using 12 lower integrals. Although very straightforward, this is numerically not 
necessarily the best method to calculate tensor integrals, since it involves the inversion of a 4 x 4- 
matrix. The method presented in [30] involves the square-root of an inverse Gram determinant of 
only 3 vectors. Also these can be chosen out of n— 1 denominators, and for high n the probability 
that a phase-space point is such that all combinations lead to small Gram determinants is rather 
low. 

For low-n integrals, i.e. for n < 4, the previous statement is obviously not true, but sev- 
eral recipes and their implementations to deal with numerical instabilities exist. Notice that, 
in renormalizable gauges, r < n, so that for low n also the cost, for example, of converting 
Passarino-Veltman functions calculated following [33] to tensor integrals like above is accept- 
able. In fact, to obtain the results in this write-up, the "Alternative Passarino-Veltman-like reduc- 
tion" from [33] was used for the 4-point integrals. It can easily be predicted when this method 
fails, in which case the method from [30] was used. For the 3-point functions, conventional 
Passarino-Veltman reduction was used. 

The end-points of the tensor reduction are scalar integrals. Also these can recursively be 
reduced further, and eventually be expressed in terms of 4-point scalar functions. For the ap- 
plication in this write-up, the unitarity-approach as presented in [11] was used to express scalar 
functions into 4-point functions directly. This choice, of course, is not in correspondence with 
the "recursivity philosophy", and in fact it strictly speaking increases the asymptotic computa- 
tional complexity to (9(n 5 2 n )@ In practice, however, it appears to be rather numerically stable 
and efficient, in particular because the formulations of [30] and [11] are compatible to large 
extend, avoiding the re-calculation of some overhead. Furthermore, for the scalar-functions, no 
numerator functions have to be evaluated, and only the coefficients for the 4-point functions have 
to be calculated, avoiding the computationally more challenging issues coming with the method 
of [1 1] in general. 

The scalar one-loop 1 -point, 2-point, 3-point and 4-point functions, finally, were evaluated 
with OneLOop [36]. Also the tensor 2-point functions were evaluated with this program. 

3 Tensor symmetrization 

Also for high-rank tensors, one has to deal with tensor-contractions in the end. Contracting rank- 
r tensors with 4 T tensor components seems hopeless, and the solution to this problem is tensor 
symmetrization. Here, we use the fact that the tensor integrals Eq.© are symmetric, and a tensor 

2 Roughly speaking C(2 n ) scalar functions need 0(n 4 ) coefficients for the 4-point functions, each of which 
involves the evaluation of 0(n) denominators. 
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integral of rank r has only 

'r + 3\ 6 + llr + 6r 2 + r 3 



/ r (6) 

r J 6 

independent components. So for symmetric tensors, the issue of contraction does not look hope- 
less at all. We only have to make sure we can calculate the symmetrized coefficients to be con- 
tracted with the tensor integrals directly. By symmetrization we mean adding tensor-components 
which are multiplied by the same tensor integral together, so 

T {1 ,2} _ T l ,2 I T 2,1 T {1 ,2,2) _ T l ,2,2 , T 2,1 ,2 , T 2,2,1 

■t=2 — x—2. ' >r=2 > 'r=3 — 'r=3 "•" ' T =3 "T 'r=3 

T {1 ,2,3} _ T 1 ,2,3 I t 2 ,3,1 I y3,l ,2 i t3,2,1 i T 2 ,1 ,3 i T 1 ,3,2 n , 

't=3 — 'r=3 "T 't=3 "T 'r=3 ' 'r=3 'r=3 ' r =3 > W 

etc.. In general, a symmetric tensor j^ 1 ^ 2 f ran k y with 4-dimensional indices can be 
represented as 

-pfvi -V2----v r } _ cr /o\ 
't — °no ,n, ,ti 2 ,n. 3 vo; 

where is the number of indices referring to dimension \i. These numbers satisfy n + ni + 
i~i2 + n.3 = r. Now suppose we have a linear recursive relation between tensors of the type 

with = K^. The solution is a product of the components of the vectors Ki to K r . To calculate 
the symmetrized product, we can cast the relation in the form 

cr cr— 1 txO _|_ rr- 1 1 

J no ,rti ,rt2 ,n3 tio— 1 ,ni ,TV2 ,ri3 r ' no ,rti — 1 ,ri2 ,TL3 r 

+ ^no.Ti! ,n2-1,n 3 + ^no.Ti! ,n2,n 3 -l ) (10) 

with the convention that n n2 n3 is identically zero whenever any of the indices is negative. 

The relation expressed by Eq.© seems rather trivial, but we will see that the tensor com- 
ponents we have to calculate satisfy very similar relations. The main difference will be that the 
simple multiplications on the r.h.s. will be replaced by more complicated contractions. This does 
not have any influence on the possibility to calculate symmetrized components directly, nor on 
the asymptotic computational complexity. In fact, like the calculation of the tensor integrals, also 
the calculation of the tensor components does not increase the asymptotic complexity given in 
Eq.©. This stems from the facts that the number of operations to be performed to calculate a 
tensor given the lower tensors is constant, and that the number of tensors entering the recursive 
equation is equal to the number of tensor integrals, i.e., no intermediate "auxiliary" tensors have 
to be calculated. The equivalents of the vectors above will essentially consist of tree-level 
off-shell currents, which are computed at a cost of 0{n 4 ). 



4 Ordered gluon one-loop amplitudes 

In the following, we wil derive recursive relations for the tensors to be contracted with the tensor 
integrals in order to arrive at one-loop amplitudes. First we repeat the known tree-level relations 
to introduce some notation. 
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4.1 Recursive relations for ordered gluon tree-level amplitudes 

The recursive relations for tree-level gluon off-shell currents are given by [40] 



j-1 j-2 j-1 

y~ ^p(Pt,k)Pk+i,j)A^ k A^ +1 ■ + y~ y~ W^ po .A^ k A^ +1 jL Af +1 j 

.k=i k=i l=k+l 



with 



and 



VU V : , p 2 ) = [ ( Pl - p 2 ) ^g„ p + ( Pl + 2p 2 ) v g£ - (p 2 + 2 Pl ) p ] 



V2 



Igs 



(11) 



(12) 



(13) 



'vpcr t "[^SpQ'Vcr 9-v9pcr g^gp-v] • 

The starting points A^ = of these recursive equations are the polarizations vectors of the 
external gluons, and we denote 

i 

Pi,j=^PTc (14) 

k=i 

where p k is the momentum of gluon k. For i < j we define p t j = 0. If pi >n = 0, then 

yt n (l,2,...,n) = g^e^p^A^ (15) 

is the tree-level color-ordered amplitude for gluon 1 to n. The full tree-level amplitude for the n 
gluons is then given by [39] 



M n = Y_ TrCT^T^ 2 ' •••T a -(-5)yL n (7t(l),7t(2) ) ...,7t(n)) , 



(16) 



where ai , a 2 , . . . , a n are the color indices of the gluons and T a are the generators of SU(N C ). 
The sum is over all permutations of the gluons except the cyclic permutations. The off-shell 
currents satisfy g^^p^- A^. = 0, and the three-point vertex can be reduced to 



V£ p (pi ,p 2 ) = ^[ (pi - p 2 rg. P + 2p 2 . g£ - 2p lp g^ ] 
The recursive equation may be represented diagrammatically by 



(17) 




(18) 



Even for a diagrammatic representation this formula is rather rudimentary, but it encodes enough 
information for our purpose. For a more detailed description of the recursive relation, we prefer 
to refer to Eq.(fTTl) instead of dressing up the diagrammatic representation. 
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4.2 Recursive relations for ordered gluon one-loop amplitudes 

A so called color decomposition as in Eq. (fT6l) also exists for one-loop amplitudes [41], and is 
given by 

M^ 1 ' - Y_ TrfT^'" •••T Q -'-')yii 1 1) (7r(l),...,7r(n)) (19) 

7TgS n /Z n 

|n/2J+1 

+ Y_ Y- T^T "' 11 ■■■T a "< c - 1 ))Tr(T Q «< c ) ■ ■ -J^) A^{n{]), . . . ,n{n)) , 

C=2 7teS n /Sn;c 

where S n;c is the subset of S n that leaves the corresponding double trace structure invariant. 
The objects An are called primitive amplitudes. They only receive contributions from diagrams 
with a particular ordering of the gluons. The partial amplitudes A n can be calculated as linear 
combinations of permutations of the primitive amplitudes. 

Given the definition of the primitive amplitudes, one can write down a recursive relation for 
off- shell currents from which the primitive amplitudes can be constructed following a relation 
like Eq.(fT5l). The blobs in the diagrammatic equation Eq.(fT8l) represent off-shell currents con- 
sisting of sums of tree-level diagrams. We represent off-shell currents consisting of diagrams 
containing exactly 1 loop by a blob with a hole, and we have 




Concerning the first line, it is clear that, since the result may only consist of one-loop diagrams, 
exactly one blob with a hole must be connected to a vertex. Since we are considering ordered 
amplitudes, all distributions of the blob with a hole over the different legs of a vertex have to be 
represented separately. 

The actual loops are generated in the second line of Eq.(l20l). Both the second and the third 
term on this line have to be added explicitly because of the ordering. These loops are constructed 
from tree-level off-shell currents with one auxiliary gluon with momentum and polarization vec- 
tor, say, q and £^(i) respectively. The index i runs from 1 to 4 such that 

4 

^^(i)e p (i) = g^. (21) 

This gluon is supposed to be virtual, and thus off-shell. The polarization vector has no real 
physical meaning, and just plays the role of the end-point of the gluonic line. We will now 
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introduce objects G^(q) including this auxiliary gluon through the formal relations 




P 2 - 



^V^(-q-pi, ir . 1 ,q+p 1J )G^(q) 

17T Z ' J 



j-1 



d 4 q 

in 2 



(22) 



k=i 



d 4 q 



> pa 



17T 



2 ^1<+1J 



The seemingly superfluous momentum shift pi,t i in the first line is to make sure that only 
inverse denominators of the form (q + p 1,1c) 2 appear in the calculation, and not for example 
(q + P2,k) 2 - This also means that the auxiliary external gluon in G* ^(q) is carrying momentum 
q + p 1,1-1 instead of q. At this point, the question is how to assign a meaning to the relations 
above, and we will explain this in the following. 

Since we are interested only in the contribution of ordered one-loop diagrams, the auxiliary 
gluon with momentum q + pi,i_i must be the first one i for every off-shell current G^(q), so 
these off- shell currents satisfy 



More explicitly, the relation is 




(23) 



G#(q) = 



,q + pi,j) 



j-i 



Y_ Vk(q+Pi,, w Pi fr Hj)G&(q)A£, 1J 



k=i-1 



j-2 j-1 

+ 11 WVG5(qK +1)l Af+i,j 

k=i-l l=k+1 



(24) 



Notice that the sum over k starts with k = i — 1 : the case that G^(q) does not contain any of 
the gluons i to j and for which it is given by 



G&-i(q) = 9 



(25) 



for every i. Introducing the symbol 



(26) 



3 Or the last one, but we choose the first one. 
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we can separate the q -dependent part of the 3-point vertex and write 

i-l 



j-2 j-1 ] 

+ 11 W^ pa G^(q)AP +1>l Af +1)j . (27) 

k=i-l l=k+l J 

From these recursive equations, we can see that G^-tq) can be expressed as follows 

G « Cq) = £ IS^(D) ^ q " -; q " (28) 

Dc{i-1,t,..,j} r=0 lljeDlM+PljJ 

where |D| is the number of elements in T>, which is a subset of the set {i — 1 , i, . . . , j} containing 
at least i — 1 and j. The tensors 9^v 2 ...-v T (D) do not depend on q. As an explicit example, we 
can write 



^ (n] _ 9^(1,4) + gU4) 
,2 ' 4lqJ " (q + Pi,i) 2 (q+PM) 2 



| 9^( 1,2,4) + 9gTI , 2, 4) q" + 9^ 2 ( 1 , 2, 4) q^ 

(q +Pi,i) 2 (q +p 1 , 2 ) 2 (q+p 1 , 4 ) 2 
9 A ^(1 ,3,4) + ffiU 3,4) q" + 9^ 2 (1 , 3,4) q"' q^ 

(q+Pi,i) 2 (q+Pi, 3 ) 2 (q+Pi,4) 2 
9^(1,2,3,4) + 9^(1,2,3,4) q-+ 9^(1,2,3,4) q^V 2 



(q +P!,i) 2 (q +P!,2) 2 (q + p ll3 ) 2 (q + pi,4) 2 

9^(1,2,3,4) q^q^q^ 

(q + pi )1 ) 2 (q +pi, 2 ) 2 (q +Pi, 3 ) 2 (q + pi, 4 ) 2 



With this observation, we can assign a meaning to the relations of Eq. (|22l) as follows. Given the 
tensor integrals 

d D q q?q?"-q? 



(30) 



i7t D / 2 n jeB [(q+pio) 2 +io] ' 

for which the numerator only contains 4-dimensional components of q, we define the object 

|D|-1 

Dc{i-1 t=0 

which does not depend on q, and assign 

5lw: pff ^A[ Hj , -f^ = -JllW!; pff A7 )k G^ 1)j . (32) 
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This will lead to a calculation of the one-loop amplitude within the scheme of [37]. We cannot 
use the tensors 9^2-%. I'D) directly to define the first line of Eq.(l22l). and we will discuss this 
below. 

First, however, we need to answer the question how to calculate the tensors S*^, ..^(D). 
Obviously, from Eq.dTTl) we can derive recursive equations for them. Writing T> = {£>', k, j}, so 
the largest two elements of D are {k, j}, we find 



i-1 

= -iS^ 2 .... r (D',k) ' 



V^ p (p 1)k)Plc+lj -)AP +1)j + Y_ W V^ +1 ,A+U 



l=k+1 



- i^ 2 ...^{D'^)X^ p A p k+ ,., (33) 

where the third line is absent for the case r = 0, and the second line is absent for the case 
r = \D\ — 1 . As the starting points of the relations we define 

S^(i) = g^ , g^(t)=0, (34) 

for any i = 0, . . . , n. 

Let us address the discussion about the computational complexity in Section [3] and compare 
Eq.(l33l) with Eq.©. The first difference is that Eq. (|33l) has tensors of rank r also on the r.h.s.. 
Secondly, Eq.(l33l involves the contraction with index "v instead of a simple multiplication. These 
differences do not influence the asymptotic computational complexity. Finally, all objects calcu- 
lated using the recursive relation are needed in Eq.(l3"TT). and no auxiliary tensors show up whose 
calculation could influence the asymptotic computational complexity. 

In order to deal with the first line of Eq. (|22)) . we introduce the objects 



Hft(q) =V^(-q-p lii _ 1 ,q+p lj -)G^(q) . (35) 
With the help of the symbol 



Y^ p = -2g^ P + g£g w + g^gpa ] (36) 



we can separate the q -dependent part of the 3-point vertex again and write 

Hft(q) = [V^(-pi,w,pij) + Y^q°] G?J(q) , (37) 

and express 

[D| 

x— x— q Vl cT 2 ■ ■ ■ q' VT 

with 

^..^(D) = V^ p (-p lil _ 1> p 1 ,,)S^ 2 ... VT (lD)+Y^ p g^... VT _ 1 (a)) . (39) 
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The second term on the r.h.s. is absent for the case r = and the first one is absent for the case 
r = \T>\. Now we simply assign 



IDI 



(40) 



Dc{i-l,i,...,j} r=0 



We can now write down the diagrammatic relation Eq.(|20l) explicitly. Denoting a one-loop 
off- shell current by B^, we have 

-' f M 

Pi -j I k=i 

j-2 j-1 



lc=i l=k+1 



J-1 



+ fi£ + £ A£ +1 j + A^ ( 



(41) 



So the program to calculate these one-loop off- shell currents is to 

1. calculate the tree-level off-shell currents A{^ with recursive equation Eq.(fTT1); 

2. calculate the tensor integrals ^pi^-^r defined in Eq.OUl): 

3. calculate the tensors 5^% 2 ...- Vr {'^ > ) using Eq.(l33l) and the tensors JC^, V2 ... Vt (2)) using Eq.(l39l); 

4. calculate the objects H^- following Eq.(l40l) and following Eq. (|3"TT) : 

5. solve Eq.([4~TT) recursively. 

The one-loop amplitude finally is given by 

= g^elv 2 n B» n _, . (42) 



4.3 Including ghost and quark loops 

In order to arrive at gauge-invariant amplitudes, ghost contributions have to be included. Also, 
one could want to include quark loops. Then, Eq.(l20l) has to be extended to 




• (43) 



11 



The factor 2 for the ghost loop is needed to arrive at gauge-invariant ordered amplitudes. The 
calculation of the ghost-tensors is rather trivial once the calculation of the gluonic ones are un- 
derstood. We prefer to focus the attention on the quark loops. The main difference is that for 
these, the q-dependence of the numerators in the off-shell currents with an auxiliary quark comes 
from the quark propagators instead of the vertices. The off-shell currents, or "off-shell spinors", 
with an auxiliary quark satisfy 




11,,'ur -<><- . (44) 



where x, y denote the spinor-indices. More explicitly, the relation can be written as 

j-i 

U*f(q) = 7— 1 —- 2 Y_ (W:(l + Pi/ U ^(q)^ +1 , j (45) 
where implicit summation also over spinor-indices is understood. Here, we introduced the matrix 

1^ = ^=7^- (46) 

The q-dependent part on the r.h.s. can be isolated even more straightforwardly than in the gluon 
case. We write 

|D|-1 

u-,?(q)= Y_ L^WP) nTq+pS (47) 

Dc{i-i,i,..,fl r=o lljeDiq+PioJ 

and find that 

U$V-vr(B) =U^ 2 .... T (23 / ,k,j) = iU^ 2 ..., r (2)',k) (r w )5p^A£ +1J 

The contribution from the quark loops to the integrands of the one-loop gluon off-shell currents, 
without the the factor —1 /N c , are given by 

Fft(q) = [q+VwV (r-v%xUjf (q) . (49) 

Remember that U*- U (q) already contains the denominator of the propagator factor on the r.h.s.. 
Now we write 

W = L Z^V-^^ n (q+P% (50) 

D C {i-1 ) i,..J}r=0 HjeDiq+PljJ 



with 



SU..^) = Pu- 1 (r^),,xU^ 2 ..., r (D',k) 

+ (r v /) y)X U^ 2 ... Vr _ 1 (D',k). (51) 
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The integrated contributions from the quark loops to the one-loop gluon off-shell currents are 
then given by 



Dc{i-1,i,...j} r=0 



(52) 



4.4 The ^ 2 -term 

Finally, to obtain gauge-invariant results, the ft 2 -term [13, 38] has to be included. The necessary 
extra vertices are given by 



K(V) = %(l/2 + A H v + N f /N c )p^g 



•vpcr 



2 n^ 

-v 



(53) 



(7/4 + A HV + 2N f /N c ) [ ( Pl - V iV^ P + 2p 2 „ g£ - 2p lp ] (54) 



-^[(6 + 2AHv + 5N f /N c )g^ 

-(5/2 + A HV + 3N f /N c ) (g^g pa + g£ 9pv ) ] , 



(55) 



where Nf is then number of quarks, and A H v = 1 in the 't Hooft-Veltman scheme and A H v = in 
the FDH scheme. The off-shell currents containing graphs with exactly 1 such vertex satisfy 
the recursive equation 



—l 



j-2 j-1 



+ y~ y~ w ^ P a R i,k A k+i ,i A r+i ,j + R k+i ,i A f+i j + A ^k A k+i ,i Rf+i j] 



k=i l=k+1 



j-1 



+ P£(pJ,)A£ + ^ V^ p (p 1)k)Pk+1J )A- k A 



p 

lc+lj 



k=i 



j-2 j-1 
k=i l=k+1 



The D^-terrn is then given by 



$-2 - 9^£nPn R ll,n-1 



(56) 



(57) 



5 Results 



The presented algorithm has been implemented in a Fort ran7 7 program which, first of all, re- 
produces all the numeric results given in [2] for multi-gluon amplitudes up to 1 gluons to at least 
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4 decimals precision (Appendix lAl). Results involving one massless quark-loop are presented in 
Appendix |Bl 

Secondly an analysis of the accuracy like in [2] and [3] has been performed, which makes 
use of the existence of a simple formula for the divergent part of color-ordered one-loop gluon 
amplitudes within dimensional regularization [42]. It is given in [2] as 



(poles) _ 



(47r) e r(1 +e)r 2 (l 



16tt 2 



HI -2e] 




5> 



-Pi,i+i 



n 

T 



A^ 



(58) 



where e = (4 — D )/2 and D is the dimension, \± is the dimensional scale. Since the divergent part 
of the one-loop amplitude should also be obtained by using the divergent parts of the initial scalar 
functions as the starting points in the calculation of the tensor integrals, there is the opportunity 
to compare the two and assess the accuracy of the latter. The left of FigQ] gives the distribution 
of the quantity 



10 



log 



A 



(1/e; 



A 



(1/e) 



A ( 2 /e] <m 

where yi! 1 1 ^ £ ^(l58l) refers to the coefficient of 1 /e in Eq.(l58l). and A^ e ^ refers to this coefficient 
calculated with the program presented in this write-up. The distribution is obtained from a large 
sample of uniformly distributed phase-space points with the same kinematical cuts as in [2, 3] 
being 

foil < 3 , p± >i >0.01 v / i , ^/ictH-cbjp + lru-ri^XM, (60) 

where r|t is the rapidity of gluon i, pj_ t t its transverse momentum and r\i its azimuthal angle, 
all with respect to the axis of the incoming gluons. Also helicity configurations where sampled 
uniformly distributed, only avoiding configurations for which the tree-level amplitude vanishes. 
The distributions for a total number of 6, 8 and 1 gluons are shown for calculations at the double 
precision level. We see a behavior compatible with [2] at the double precision level, and slightly 
better than [3] at double precision level. The lower graph at the left of Fig{T] shows the same 
distributions in a log scale for the y-axis in order to highlight the right tail. 

Obviously, for a small part of the phase space points the accuracy becomes unacceptably bad. 
Results of re-evaluation for these at higher precision are given on the left of Figf2l Presented are 
the right tail of the distribution for n = 10 starting from —4, and distributions obtained with the 
same set of phase space points for different options of applying quadruple precision arithmetic. 
Evaluating the tensor integrals only (but including the scalar functions) at quadruple precision, 
the tail is shifted to the left to a large extend, but still phase space points may show up leading 
to an unacceptably bad accuracy. The more expensive option of evaluating everything at the 
quadruple precision level, however, moves the whole tail below —1 6. 

Another hint at the accuracy of the program can be given by the extend at which gauge 
invariance is satisfied. A one-loop amplitude should vanish whenever the polarization vector of 
any of the external gluons is replaced by the momentum of that gluon. The quantity 



(59) 



10 



log 



Re^ln (£t <- Pi 



+ 



lmA^ ] 



M 



lmA\ 



(61) 
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-16-14-12-10 -8 -6 -4 -2 2 -16-14-12-10 -8 -6 -4 -2 2 

1/epsilon pole gauge invariance 

Figure 1: The distribution of the quantity in Eq.(l59l) (left) and Eq.doTT) (right) for calculations at 
the double precision level. The lower graphs represent the same distributions as the upper ones, 
but with a logarithmic scale for the y-axis. 

may serve as a measure of the number of decimals being eliminated by performing such a re- 
placement in a numerical calculation, which then again may give an estimate of the accuracy. 
The right of Fig \\\ presents the distribution of (the finite part of) this quantity, obtained from the 
same sample of phase-space points as before, now however including helicity configurations for 
which the tree-level amplitude vanishes. The distributions are compatible with the ones on the 
left. The lower graph at the rights show the same distributions again in a log scale for the y-axis. 
The right tails behave similarly to the ones left, and the same holds for the distributions in Figj2l 
In Table Q] the typical cpu-times t loop are given needed for 1 evaluation of the one-loop am- 
plitude for a number of gluons from 4 to 10 on a 2.80GHz Intel Xeon processor. They are 
determined by taking the average over the evaluations for a large number of different phase- 
space points. The numbers are roughly comparable with those in [2, 3], but do show a worse 
behavior as function of the number of external particles. The numbers for t\ oov /{x\. 4 2 n ) seem 
to converge, supporting the statements in Section [3] about the computational complexity being 
0{n 4 2 n ). The numbers for ti oop /t tree , where t tree is the cpu-time needed for 1 evaluation of the 
tree-level amplitude, give a machine-independent measure of the computational cost. Notice that 
they increase as 2 n , consistent with the computational complexities of 0{n 4 2 n ) for the one-loop 
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0.1 



0.01 



0.001 




0.1 



0.01 



0.001 



-20 -15 -10 -5 
1/epsilon pole 




-20 -15 -10 -5 
gauge invariance 

Figure 2: The right tail of the distributions for n = 10 of Fig |T| (evaluated at double preci- 
sion level), and the distributions for the same phase space points when the tensor integrals are 
evaluated at quadruple precision, and when everything is evaluated at quadruple precision. 



n 


4 


5 


6 


7 


8 


9 


10 


tioop(ms) 


2.762 


10.15 


34.37 


109.8 


335.1 


965.2 


2744 


t loop (p.s)/(n 4 2-) 


0.6744 


0.5077 


0.4144 


0.3573 


0.3196 


0.2873 


0.2680 


tloop/ftree/1 3 


0.2990 


0.6102 


1.180 


2.244 


4.104 


7.919 


15.17 



Table 1 : Typical cpu-times needed for 1 evaluation of the one-loop amplitude on a 2.80GHz Intel 
Xeon processor. 



amplitude, and 0{n 4 ) for the tree-level amplitude. 



6 Conclusion 

An algorithm was presented to calculate multi-gluon one-loop amplitudes using tensor integrals, 
which was shown to be competitive with existing programs using the unitarity-approach up to 
a number of 10 gluons. It uses universal recursive relations for tensor integrals, independent of 
the amplitude being calculated. It also uses recursive relations for ordered gluon amplitudes, 
however in this respect it can straightforwardly be generalized to any field theory by extending 
known recursive relations at tree-level [43-45] to the one-loop level, as was done for the ordered 
gluon amplitudes in this write-up. 
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A Reproduction of existing results 

This appendix contains the reproduction of the explicit numeric results in [2] up to a number of 
1 gluons. The phase-space points at which the amplitudes were evaluated can be found in [2]. 
Presented is the absolute value of the coefficients of e~ 2 , e~ ] and e° of the amplitudes. Results 
labelled with "tree" give the latter value for the tree-level amplitude. The upper-left box in each 
table gives the helicity configuration. The results from [2] are labelled with " [2]". 

A.l n = 6 



++++++ 


e- 1 


e-' 


e° 


tree 
tree [2] 
unit [2] 
num [2] 

tnsr 


0.0000000000000000E+00 
0.1060660419488780E-13 
0.1349521124077591E-10 


0.0000000000000000E+00 
0.3813284749527035E-13 
0.779228371397191 IE- 10 


0.2215923815877299E-14 
0.1767767365814634E-14 
0.5298064836438550E+00 
0.5298064836612950E+00 
0.5298064837199270E+00 




-+++++ 


e- 2 


e-' 


e° 


tree 
tree [2] 
unit [2] 
num [2] 

tnsr 


0.101 125576124171 1E-10 
0.2377895374324842E-12 
0.6401972138850929E-10 


0.6753625348984687E-09 
0.8549005883762705E-12 
0.1571282783053386E-08 


0. 1 947 1 241 36075 826E- 13 
0.3963158957208070E-13 
0.3259967043518990E+01 
0.3259967054272360E+01 
0.3259967053973336E+01 




- - + + ++ 


e- 2 


e-' 


e° 


tree 
tree [2] 
unit [2] 
num [2] 

tnsr 


0. 1709476899026590E+03 
0.1709476899026590E+03 
0.1 709476899026 128E+03 


0.6145908783763959E+03 
0.6145908783763970E+03 
0.6 14590878375069 1E+03 


0.2849 1 28 1 65044324E+02 
0.2849 1 28 1 65044320E+02 
0. 1 373747535008540E+04 
0. 1 373747535008280E+04 
0.13737475350071 19E+04 




— 1 h -+ 




e- 1 


e° 


tree 
tree [2] 
unit [2] 
num [2] 

tnsr 


0. 1 883229237004670E+02 
0.1883229237004850E+02 
0.1883229236988423E+02 


0.6770582934748300E+02 
0.6770582928695769E+02 
0.67705 8292866730 1E+02 


0.3138715395008085E+01 
O.3138715395008O8OE+01 
0.1510439503289600E+03 
0.1510439503379470E+03 
0.1510439503340157E+03 
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H 1 h- 


e- 1 


e- 1 


e° 


tree 
tree [2] 
unit [2] 
num [2] 

tnsr 


0. 1 883229237005540E+02 
0. 1883229237004850E+02 
0. 1883229236987900E+02 


0.6770582928570479E+02 
0.6770582928695769E+02 
0.67705 829287 1 7265E+02 


0.3138715395008091E+01 
0.3138715395008080E+01 
0.1537801015298360E+03 
0.1537801014159860E+03 
0.1537801016156306E+03 



A.2 n = 7 



+++++++ 


e- 1 


e- 1 


e o 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.0000000000000000E+00 
0.0000000000000000E+00 
0.6861582055762283E-11 


0. 1 256534542409480E-09 
0. 12501701 11559883E-14 
0.3862018590575813E-10 


0.2703922191931054E-14 
0.0000000000000000E+00 
0.3101695329720260E+00 
0.3101695334831830E+00 
0.3101695335777539E+00 




-++++++ 


e~ 2 


e- 1 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.36782128743 19657E- 12 
0.2713533399763100E-14 
0.4384180396142481E-10 


0.7209572152581734E-12 
0.8924875 144594874E- 14 
0.2104965953958677E-09 


0.1648597081617964E-14 
0.0000000000000000E+00 
0.1920528148108100E+00 
0.1920528147653950E+00 
0.1920528150979991E+00 




h + + + + 




e- 1 


e o 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0. 14746289842 1 6 140E+02 
0. 14746289842 1 6 140E+02 
0.1474628984215924E+02 


0.4850089396312140E+02 
0.4850089396312130E+02 
0.4850089396308871E+02 


0.2106612834594487E+01 
0.2106612834594490E+01 
0.873 1521551387900E+02 
0.873 1521551386510E+02 
0.873 1521551379141E+02 



-+-+-+- 


e~ 2 


e- 1 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.7713059766610930E+00 
0.7713059766610950E+00 
0.7713059766434806E+00 


0.2536843489960730E+01 
0.2536843489960750E+01 
0.2536843489868709E+01 


0.1101865680944418E+00 
0.1 101865680944420E+00 
0.5933610502945470E+01 

0.5933610502629016E+01 




H h-H h 


e- 2 


e-' 




tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.77 1 30597666 1 0930E+00 
0.77 1 30597666 1 0950E+00 
0.7713059766436269E+00 


0.2536843489960740E+01 
0.2536843489960750E+01 
0.2536843489862686E+01 


0.1101865680944424E+00 
0.1101865680944420E+00 
0.6042012409916140E+01 

0.6042012409615074E+01 



A.3 n = 8 



++++++++ 


e~ 2 


e-' 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.0000000000000000E+00 
0.3853462894343397E-14 
0.1336896918161384E-10 


0.0000000000000000E+00 
0.1441159379540454E-13 
0.565240309721661 IE- 10 


0.2961844168510185E-15 
0.0000000000000000E+00 
0. 1 9670060069569 1 OE+00 
0.1967006007382010E+00 
0. 19670060 1 2222209E+00 




-+++++++ 


e- 2 


e- 1 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.0000000000000000E+00 
0.1805821909003967E-13 
0.2640076774290492E- 10 


0. 1 965638 1 04048654E-09 
0.6753606439965886E-13 
0.1008615996797149E-09 


0.295 1926495253574E- 14 
0.2257277386254959E-14 
0.5287747 164930630E+00 
0.5287747176521700E+00 
0.5287747 173953582E+00 
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■ — h + + + ++ 


e- 1 


e- 1 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.3466551355735610E+02 
0.3466551355735680E+02 
0.346655 1 355735630E+02 


0.1296458053471450E+03 
0.1296458052914090E+03 
0. 12964580529 1 3800E+03 


0.43331891 946695 86E+01 
0.4333189194669600E+01 
0.2742997734349260E+03 
0.2742997734349000E+03 
0.2742997734347 1 39E+03 




- + - + -H h 


e- 2 


e-' 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.5809218091 107730E+00 
0.5809218091 108460E+00 
0.5809218091244186E+00 


0.2172593680235970E+01 
0.2172593682447690E+01 
0.2172593682447645E+01 


0.7261522613885360E-01 
0.7261522613885579E-01 
0.54763038 19766790E+O1 

0.5476303819972849E+01 



H 1 1 h- 


e- 2 


e-1 


e o 


tree 






0.7261522613885418E-01 


tree [2] 






0.7261522613885579E-01 


unit [2] 


0.5809218091 108620E+00 


0.2172593687810420E+01 


0.4925500546307290E+01 


anly [2] 


0.5809218091 108460E+00 


0.2 172593682447690E+0 1 




tnsr 


0.5809218091 153071E+00 


0.2172593682480466E+01 


0.4925500546470287E+01 



A.4 n = 9 



+++++++++ 


e~ 2 


e" 1 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.48602698362923 1 6E- 1 1 
0.2693624076029 1 1 6E- 1 2 
0.1242803468718531E-06 


0. 1 845 1 93695700690E-07 
0.1176695244346755E-11 
0.1395624387049265E-05 


0.3090869336705567E-13 
0.2992915640032351E-13 
0.5666555617062950E+01 
0.5666555580473110E+01 
0.5666561857344509E+01 




-++++++++ 


e- 2 


e- 1 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.3938371378126140E-10 
0.8202679137223861E-12 
0.1761715389271764E-05 


0.2340429860576292E-07 
0.3583296428617654E-11 
0.2085999886200371E-04 


0.7191085595712448E-13 
0.9114087930248735E-13 
0. 1 0620864606 14280E+0 1 
0.1062086467981750E+01 
0.1062166559796595E+01 




+ + + + + + 


e- 2 


e-' 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.2909067010969220E+03 
0.290906701 1509570E+03 
0.2909066975088621E+03 


0.1270810334861320E+04 
0.1270810336301850E+04 
0.1 27081 0392 170893E+04 


0.3232296679455 183E+02 
0.3232296679455080E+02 
0.3625430616705210E+04 
0.3625430616705940E+04 
0.3625430209440870E+04 




— 1 1 1 h - 


e- 2 




e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.408 1 69769666 1 860E+01 
0.408169769731 1550E+O1 
0.4081697504765963E+01 


0. 1 783067767208 1 40E+02 
0. 1 783067768078440E+02 
0.1783065940829036E+02 


0.4535219663678330E+00 
0.4535219663679500E+00 
0.57 10639504628740E+02 

0.5710625504238800E+02 




H I-- + -H h 


e- 2 


e-' 


e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.4081697696620550E+01 
0.408169769731 1550E+01 
0.408169755840971 1E+01 


0. 1 783067764548420E+02 
0.1783067768078440E+02 
0.1783066268732069E+02 


0.4535219663678577E+00 
0.4535219663679500E+00 
0.5501538077075760E+02 

0.550153793 1596933E+02 
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A.5 n = 1 



++++++++++ 


e- 1 


e- 1 




tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.2616999209810146E-12 
0.7645214091184737E-12 
0. 1 199729722684045E-06 


0.7453 142378465002E-06 
0.3853184186191476E-11 
0.1810959673219180E-05 


0.8422572777655544E-13 
0.76452 1 409 1184737E- 13 
0.18434901 12846700E+02 
0. 1 843490 1 1 28467 1 OE+02 
0. 1 843487909054984E+02 




-+++++++++ 


e- 1 




e° 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.1729567134060808E-10 
0.3138928592085274E-11 
0.9039805077823879E-06 


0.3462486730362966E-05 
0.1582018484813023E-10 
0.8416576413616293E-05 


0.15381906621 18770E-12 
0.3138928592085274E-12 
0.1411 806902836740E+02 
0.1411 806902836920E+02 
0.141 1799961 139769E+02 



+ + + + + 


e- 1 


e-' 


e° 


tree 






0.4899726956663458E+03 


tree [2] 






0.48997269566634 1 0E+03 


unit [2] 


0.4899726956656070E+04 


0.2469460004000990E+05 


0.7584491014580890E+05 


anly [2] 


0.48997269566634 1 OE+04 


0.2469460004768270E+05 


0.7584491014578140E+05 


tnsr 


0.4899726958888808E+04 


0.2469460005 194075E+05 


0.7584491017379209E+05 



- + - + -H 1 h 


e- 1 


e- 1 


e" 


tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.93461 13719987591E+02 
0.93461 13720089021E+02 
0.9346 1 1 3676 1 27024E+02 


0.47 1 0436787027 1 10E+03 
0.47 1 0436772479390E+03 
0.4710436797756468E+03 


0.93461 13720088734E+01 
0.93461 13720089020E+01 
0. 1 48 1 274476056640E+04 

0.1481 274475543962E+04 




H 1 1 1 h- 




e-' 




tree 
tree [2] 
unit [2] 
anly [2] 

tnsr 


0.9346113719956180E+02 
0.93461 13720089021E+02 
0.93461 13696575948E+02 


0.47 1 0436740057420E+03 
0.47 1 0436772479390E+03 
0.47 1 0436794794809E+03 


0.93461 13720088464E+01 
0.93461 13720089020E+01 
0. 1503970258031 1 1 OE+04 

0. 1 503970260993868E+04 



B Quark loops 

This appendix contains results with one massless quark-loop exclusively without the gluoninc 
loops, and the full one-loop amplitude including one massless quark-loop. The former are la- 
belled with "excl", and the latter with "incl". The phase-space points are the same as in Ap- 
pendix [A] and are given in [2]. Presented is the absolute value of the coefficients of e~ 2 , e _1 
and e° of the amplitudes. Results labelled with "tree" give the latter value for the tree-level 
amplitude. The upper-left box in each table gives the helicity configuration. 

B.l n = 6 



++++++ 


e- 2 


e" 1 


e° 


tree 
excl 
incl 


0.3444649545719067E-1 1 
0.548408240608 1639E-11 


0.3563314196395845E-11 
0.5497590035137738E-10 


0.2215923815877299E-14 
0. 1 76602 1 6 1 24 1 2397E+00 
0.3532043223779056E+00 
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-+++++ 


e~ 2 


e-' 


e° 


tree 
excl 
incl 


0.3588818407048125E-10 
0.8158777253122160E-10 


0.2958732120038221E-09 
0.1238241794955264E-08 


0. 1 947 1 241 36075 826E- 13 
0.1086655684552452E+01 
0.2173311369421295E+01 




- -++++ 


e- 2 


e- 1 


e° 


tree 
excl 
incl 


0.3083625492388671E-10 
0.1709476899026306E+03 


0.633139592231228OE+01 
0.6094559286552980E+03 


0.2849 1 28 1 65044324E+02 
0.2228503047232156E+02 
0.1363324204445743E+04 




- + -H h 


e~ 2 


e-' 


e° 


tree 
excl 
incl 


0. 828293 1375686374E- 10 
0.1883229236997337E+02 


0.6974923101314857E+00 
0.67 140 1 4235334223E+02 


0.3138715395008085E+01 
0.3642 1 289320228 1 5E+0 1 
0. 1496072786704808E+03 




H 1 h- 


e- 2 


e-' 


e° 


tree 
excl 
incl 


0.1130982098649374E-10 
0.1883229236989885E+02 


0.697492309824601 8E+00 
0.67 140 1 42354 1 9050E+02 


0.3138715395008091E+01 
0. 1 3999555504575 1 3E+0 1 
0.1526894583161683E+03 



B.2 n = 7 



+++++++ 


e~ 2 


e- 1 




tree 
excl 
incl 


0.27343548 1 8049004E- 12 
0.5758778544760509E-12 


0.1535659632517864E-11 
0.3391282384997409E-11 


0.270392219193 1054E-14 
0.1033898444969325E+00 
0.2067796889938174E+00 




-++++++ 


e- 2 


e-' 


e° 


tree 
excl 
incl 


0.7829286354785003E-11 
0.1149200607547395E-10 


0.40383785465 19624E- 10 
0.6152252106547420E-10 


0.1648597081617964E-14 
0.6401760500069484E-01 
0.128O352O99351288E+00 




H + + + + 


e- 2 


e" 1 


e° 


tree 
excl 
incl 


0.3039208000395293E-12 
0. 14746289842 1 5543E+02 


0.468 1 36 1 854653245E+00 
0.4815969760605956E+02 


0.2106612834594487E+01 
0.7724537168857577E+00 
0.8723793738684391E+02 




-+-+-+- 


e- 2 


e-' 


e° 


tree 
excl 
incl 


0.4415997997437337E-11 
0.7713059766680262E+00 


0.2448590404383259E-01 
0.2518997184786409E+01 


0.1101865680944418E+00 
0.8801500600694001E-01 
0.5849371002103936E+01 




H h-H h 


e- 2 


e-' 


e° 


tree 
excl 
incl 


0.4278136193897277E-11 
0.7713059766681208E+00 


0.2448590401493046E-01 
0.251 8997 184763565E+01 


0.1101865680944424E+00 
0.9578987033471566E-01 
0.5952686996619508E+01 



B.3 n = 8 



++++++++ 


e~ 2 


e- 1 


e 6 


tree 
excl 
incl 


0.6535721923323608E-12 
0.6936524320241993E-11 


0.2946477554703654E-10 
0.2689757765745469E-10 


0.2961844168510185E-15 
0.6556686705957517E-01 
0.131 1337340887097E+00 




-+++++++ 


e 


e" 1 


e° 


tree 
excl 
incl 


0.29935 17822283384E- 11 
0.1682824582337319E-10 


0.4138014710832325E-10 
0.2282288551418147E-10 


0.2951926495253574E-14 
0.1762582391337885E+00 
0.3525164784441872E+00 
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■ — h + + + ++ 


e- 1 


e- 1 


e° 


tree 
excl 
incl 


0.2048884758899206E-11 
0.346655 1355734642E+02 


0.9629309321605227E+00 
0.1288994364902356E+03 


0.43331891 946695 86E+01 
O.1952623627599051E+O1 
0.2738716469885952E+03 




- + - + -H h 


e- 2 


e- 1 


e° 


tree 
excl 
incl 


0.6007828400852260E-11 
0.5809218091 173581E+00 


0.1613671694024866E-01 
0.216008609576755 1E+01 


0.7261522613885360E-01 
0.7432198837943639E-01 
0.5416218143354230E+01 




H 1 1 h- 


e- 2 


e- 1 


e° 


tree 
excl 
incl 


0.2176051463652288E-11 
0.5809218091 128370E+00 


0.1613671696641898E-01 
0.2160086095780097E+01 


0.7261522613885418E-01 
0.38874483458 10546E-01 
0.4897808762927628E+01 



n = 9 



+++++++++ 


e- 1 


e-' 


e° 


tree 
excl 
incl 


0.1679041 191381272E-07 
0.9845976672993657E-07 


0.4256557696335983E-06 
0.1169508469550796E-05 


0.3090869336705567E-13 
0.1888851 178320152E+01 
0.3777703015967018E+01 



-++++++++ 


e~ 2 


e- 1 


e 6 


tree 
excl 
incl 


0.7127945940983221E-06 
0.1 1241217268071 17E-05 


0.7 14254641 8827421E-05 
0.1415224198572896E-04 


0.7191085595712448E-13 
0.3540723883 196215E+00 
0.7080938879221086E+00 



h + + + + + + 


e- 1 


e"' 


e° 


tree 
excl 
incl 


0.36971 11686288034E-05 
0.2909066988869654E+03 


0.71 828733 16896260E+01 
0.1264862610175780E+04 


0.3232296679455 1 83E+02 
0.8760538757567446E+01 
0.3621 1890401 16724E+04 



— 1 1 1 h - 


e~ 2 


e-' 




tree 
excl 
incl 


0.2296878137953301E-06 
0.408 169755745 1 384E+01 


0. 1 0077477 1 53 1 8806E+00 
0. 1 77472 1 263709833E+02 


0.4535219663678330E+00 
0.347220 1 5630873 1 4E+00 
0.5675986177672981E+02 



H I-- + -H h 


e- 2 


e" 1 




tree 
excl 
incl 


0.2530328589867O88E-06 
0.4081697604003098E+01 


0.1 007792647 137644E+00 
0.1 774721 6487 17686E+02 


0.4535219663678577E+00 
0.3151037218153409E+00 
0.54701 12910999667E+02 



n = 10 



++++++++++ 


e- 1 


e- 1 


e 6 


tree 
excl 
incl 


0.77425 16776772886E-07 
0.183364691 17 15754E-06 


0.8116698643478135E-06 
0.2136772265813164E-05 


0.8422572777655544E-13 
0.6144957172059176E+01 
0.1228991 105613530E+02 




-+++++++++ 


e- 2 


e- 1 


e° 


tree 
excl 
incl 


0.5628204459080735E-06 
0.1437262034778626E-05 


0.3925867464691803E-05 
0.8702188415939938E-05 


0.15381906621 18770E-12 
0.4705985971256853E+01 
0.9412025021303032E+01 




- - + + + + + + ++ 


e 


e"' 


e° 


tree 
excl 
incl 


0.1254733694226626E-05 
0.4899726955335932E+04 


0.1088828222285149E+03 
0.2460028096277 1 22E+05 


0.4899726956663458E+03 
0.2631225655494603E+03 
0.757 14946849044 1 6E+05 
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- + - + -H 1 h 


e~ 2 


e" 1 


e° 


tree 
excl 
incl 


0.4577998694131886E-06 
0.93461 13693450476E+02 


0.2076909542983419E+01 
0.4692445627736100E+03 


0.93461 13720088734E+01 
0.7503 126329349836E+01 
0.1474914895661399E+04 




H 1 1 1 h- 


e- 2 


e-' 


e° 


tree 
excl 
incl 


0.4345102202108285E-06 
0.93461 1370501 1550E+02 


0.2076912724815172E+01 
0.46924456 1 67038 1 9E+03 


0.93461 13720088464E+01 
0.6824080572856147E+01 
0. 1498766884054935E+04 
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